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Abstract. "We find a class of exact solutions of differentially rotating dust 
in the framework of General Relativity. There exist asymptotically flat 
space-times of the flow with positive mass function that for radii suf- 
^^ I ficiently large is monotonic and tends to zero at infinity. Some of the 

C ' ' !3 ' space-times may have non-vanishing total angular momentum. The flow 

is essentially different from another exactly solvable flow described by van 
Stockum line element. 



1. Introduction 

The paper is devoted to investigation of solutions of self-gravitating cylindrically 
symmetric dust flow on circular orbits along field lines of locally non-rotating ob- 
servers in the framework of General Relativity. The fiow is non-rigid (differential) 
and non-expanding. It turns out the resulting space-times are globally regular, 
apart from the points where singularities may be located, and that the congruence 
of locally non-rotating observers used to define energy-momentum exist globally, 
as well. Consequently, the problem of the flow is well posed. We shall call such 
a fiow the K. flow for brevity. The main drawback of the fiow is that the proper 
energy density is negative definite and that the fiow is purely relativistic - has no 
newtonian limit. 

The result of our paper is the construction of the multipole sequence of asymp- 
totically flat, stationary, and cylindrically symmetric external solutions to which a 
broad class of other asymptotically fiat solutions can be decomposed. We construct 
also the corresponding multipole sequence of internal space-times. The multipole 
representation is enumerable by a discrete parameter. There exist also a contin- 
uum spectrum of basic solutions of which we do not examine here as they are 
nonanalytic on the axis of rotation. The stationary part of multipolar solutions 
that are asymptotically fiat we call external multipoles, and the stationary part 
of multipolar solutions that are not asymptotically flat we call internal multipoles. 
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We showed that external multipoles, and any asymptotically flat space-times of /C 
flow, are massless no matter what is pasted smoothly as an alternative internal 
completion of the space-times. Despite the fact, some of the space-times can have 
non- vanishing total angular momentum. 

As in standard gravity a broad class of solutions can be decomposed to the multi- 
poles. Even though the resulting space-times are solutions of Einstein equations of 
a differentially rotating dust, thus physically more interesting than van Stockum- 
Bonnor flows PP, [3|, 0], it seems that they are of no physical concern as local 
energy density is not positive definite. Anyway, we decided to publish the solutions 
as they present another class of exactly solvable dust flows in general relativity. 
Similarly, as the exactly solvable van Stockum-Bonnor class of asymptotically flat 
flows along time Killing vector fleld, the asymptotically flat external multipoles of 
/C flows, which must have vanishing total mass irrespectively of internal deforma- 
tions^, have no newtonian limit and are useless in astrophysical situations. As far 
as exotic forms of matter concern, the internal solutions of both the models can 
be still considered as physically interesting and worth to be examined if completed 
smoothly by, say, asymptotically flat external (Kerr) vacuum solutions. 

It turns out that the structure function K = ^t] oi the flow is the same as for the 
van-Stockum Bonnor flow 0j, nonetheless, both the flows and geometry of the re- 
sulting space-times are qualitatively different. An asymptotically flat van Stockum- 
Bonnor flow is rigid, non-expanding, and has positive energy density proportional 
to the square of the vorticity scalar which does not vanish. An asymptotically flat 
/C flow is differential, non-expanding, locally non-rotating, and has negative definite 
energy density proportional to the square of the non-vanishing shear tensor. Dust 
of van Stockum-Bonnor does not rotate with respect to asymptotic observers, it has 
a point-dependent physical velocity as measured with respect to the locally dragged 
inertial frames and which is proportional to K, whereas dust of a /C flow rotates 
differentially with respect to asymptotic observers with angular velocity propor- 
tional to K, and has vanishing physical velocity. The peculiarities are possible due 
to purely relativistic effects and show the flows have no newtonian limits. 

As an example we constructed a z-symmetric, cylindrically symmetric, and sta- 
tionary for radii sufficiently large, asymptotically flat space-time which is globally 
smooth (with the exception of point-like singularities of curvature residing on the 
axis of rotation), and which can be decomposed globally to the external multipoles. 
The solution has nonzero total angular momentum and zero total mass. 

2. Setup 

We assume the global existence of the cylindrical symmetry space-like Killing 
vector T] with closed field lines, together with the global existence of time transla- 
tion Killing vector ^ of which field lines are opened, time-like for radii sufficiently 
large, and asymptotically normalizable to unity. Having in mind stationary and 
asymptotically flat space-times with cylindrical symmetry, we shall proceed as in 
the standard theory of rotating stars in the relativistic astrophysics 0. The as- 
sumption of asymptotical flatness allows for the unique determination of Killing 
vectors ^ and t] by the defining properties. In addition one assumes that asymptot- 
ically ^T] — > 0. Field lines of the Killing vectors, which clearly are frame independent 



one can imagine such a deformation as a replacement of the interior by a region of arbitrary 
space-time such the resulting space-time be still smooth 



objects, may be viewed as two of four coordinate lines in some particular coordinates 
in which the time coordinate t runs along open lines of ^ and the cyclic coordinate 
(j) along closed lines of rj, that is, ^ = dt and rj = c?^ by definition. The other 
two, denoted by p and z, are arbitrary internal coordinates in a two dimensional 
subspace orthogonal to ^ and rj. In this coordinates the most general line element 
of a cylindrically symmetric and stationary space-time is fully determined by four 
structure functions v(p,z), ip{p,z), oj{p,z) and p{p, z) and reads 2^ 

(2.1) ds^ = e^^dt^ - e"^ (d</> - uAtf - e^^ (dp^ + Az^) . 

Bardeen coordinates are distinguished by the property that in these coordinates 
Killing vectors ^ and rj attain particularly simple form 

In asymptotically flat space-time one can introduce cylindrical coordinate system 
in which the line element at infinity reduces to 

ds^ -^ At^ ~ dp2 - p2 J02 _ j^2 

Axis of the coordinates frame are attached to 'fixed stars', otherwise the axis would 
rotate. This in turn would be in contradiction with the assumption that dt can be 
asymptotically normalized to unity. One assumes, therefore, that asymptotically 
condition e^^ —up'e^'^ > should hold. The latter is not the case e.g. in a rotating 
frame of reference that rotates with constant angular velocity. Put differently, 
Bardeen coordinates are asymptotically inertial (anyway, this metric can be treated 
more formally and be used to describe solutions that are not asymptotically flat). 
As it turns out later, IC flow can be defined globally, that is, locally non-rotating 
observers exist for all points (apart from singularities), but, unlike for van-Stockum 
Bonnor flow, space-time of K. flow is not globally stationary. The K. flow is differ- 
ential (the shear tensor does not vanish). The angular velocity of asymptotically 
flat IC flow is, by construction, identical to the angular velocity of dragging of in- 
ertial frames Q. As so, dust rotates on circular orbits with respect to asymptotical 
stationary observers even though it locally does not rotate (the vorticity vector of 
the flow vanishes identically) and its physical velocity, measured with respect to 
dragged local inertial frames, is also identically zero. 

2.1. The line element of K. flow. We derive the line element of the IC flow step 
by step from the general form of the line element given by 12.11 We specify the 
structure functions uniquely by the requirement that dust space-time trajectories 
are identical with field lines of four-velocity field of locally non-rotating observers 
and that Einstein equations are satisfied. 

2.1.1. A note on locally non-rotating observers. It should be clear that irrespec- 
tively of any particular reference frame, a generic cylindrically symmetric and sta- 
tionary flow is fully determined by a four-velocity field 

u{n) = z{n){^ + nr]), 

provided 

Z{n)-^ =^^ + 2n^rj + n^rjrj > 0, £^n = = £r,n. 



In an asymptotically flat space-time fl has the interpretation of the angular velocity 
measured with respect to 'fixed stars'. In particular, the condition U{Vt)ri = gives 

o ^'^ ~ 

\l = — OJ. 

VV 

The vorticity tensor^ of velocity field n = U{Cj) vanishes identically, which can 
be checked by a direct calculation in Bardeen coordinates 12. II Therefore, the ob- 
servers are called locally non-rotating, although they move differentially on circular 
orbits with respect to 'fixed stars' with angular velocity Cj of dragging of iner- 
tial frames. Indeed, congruence of locally non-rotating observers distorts without 
changing proper volume^ 

n^(n) = 0, 0(n)EEO, a' {n) = ---'"'" 



2.„^_ 1(W) (VcZ-)^ . ill 



iv ii 
vv iv 



vv 



Local inertial frames tangent to n play the role of local standards of rest with 
respect to which 'physical' velocity is measured. Therefore, the physical velocity 
of a particle moving with four- velocity U{fl) is the relative velocity of two inertial 
observers momentarily co-moving with U{Q) and n, respectively. After Special 
Relativity we define the relative velocity as the hyperbolic tangent of the hyperbolic 
angle x between four-velocities [/(fi) and n. In Bardeen coordinates it reads 

(2.2) tanhx^e'^-^in-Cj). 

The result agrees with that of |5] where more 'physical' derivation was presented. 

2.1.2. Specifying the line element. The energy-momentum tensor of dust matter 
moving along field lines of vector field n reads T — T>U ® U where 

[/((I-) = Z(J)) (I + WT7) ^ C/^ = e"^[l, 0,^,0]. 

From the previous construction it follows that dust in asymptotically flat JC flow 
moves differentially on circular orbits with angular velocity of dragging of inertial 
frames uj as seen by an asymptotic stationary observer. This flow has vanishing 
vorticity vector, that is, it does not rotate locally, and has no physical velocity 
- it is at rest with respect to congruence of local standards of rest dragged with 
angular velocity u). It follows also that angular momentum per mass element of the 
flow —If^T]^ is identically zero as [/ = n and ry are orthogonal. Nevertheless, as 
we shall show later, total angular momentum does not vanish for particular kinds 
of asymptotically flat space-times of JC flow which must contain singularities that 
contribute total angular momentum of the space-times. 

Einstein's equations G = kT and the contracted Bianchi identities V^G^jy ~ 
yield local conservation law V^T^^ = that for dust gives continuous flow along 



see footnote U] 

Dilation tensor &{u) and (traceless) shear tensor cr{u) are defined for a velocity field u 
as Qfii, = S7c,u°'hfiv and o-pi^ = V(a«^)/i",,/i ^, — \SJaU°'hp,u, where h'^^ = S^^ — ui^Uu is 
projector h{u) onto subspace orthogonal to u. Vorticity tensor f2(tt) for u is defined as u)^^ = 
^ [a'>^ l3]t'-°' u,^ ,j ^i^d it yields a derivative quantity that characterizes vertex - the vorticity vector 
ujt^ = ^ "^ —— UdUi^j. One defines also square of dilation scalar cr'^{u) = ^cr^^cr^i^ > and square 
of vorticity scalar fi-^(ti) = itj^'aj^jy > 0, then fl^{u) = — Li;''ajp. 



geodesic paths 

V^ (VW) = and [/"V^f/^ = 0. 

The continuity equation is satisfied identically for /C flow on the power of Killing 
equations and of symmetries of energy density ^V = and r/D = 0. In Bardeen 
coordinates IfV^Uf^ ~ {0, —dpX, 0, —dzX} which implies that A must be a constant. 
Another constraint on structure functions follows from the linear combination of 
Einstein's equations 

The left side is identically zero on the power of Einstein equations, while the ex- 
pression to the right holds for pressureless perfect fluid in Bardeen coordinates. We 
have not yet specified the arbitrary coordinates p and z. Asymptotical flatness in 
cylindrical coordinates requires that e^^ — > 1 and e'^'^ -^ p^, thus e^^'^' -^ p. The 
simplest way to concord the facts is to choose as a solution of the above equation 

(2.3) A = 0, e"^ = p. 

We have thus shown that /C flow requires at most two structure functions w and p. 

Killing vector rj is always space-like rfT] = —p^ thus the space-time of IC flow is 
globally cylindrically symmetric. As ^^ = 1 — w^p^ it is also stationary in regions 
where uj'^p'^ < 1. Note that det[5] = — p^e'*'^ does not depend explicitly on ui thus 
the boundary surface w^p^ = 1, which demarcates stationary and non-stationary 
regions of a space-time, is nonsingular, and in a sense, is alike to Schwarzschild 
horizon. 

We stress the construction of the energy-momentum tensor of IC flow makes sense 
globally as C/^f/^ > 0, or ^^ + 2a)|r; + LJ^r/q = 1 > for /C flow. C7(a)) is time-like 
even though in some regions the resulting space-time may be non-stationary, and 
so, C/(w) be a combination of two space-like vectors. It follows that in space-times 
of K. flow the locally non-rotating observers exist globally. 

Similarily as in jjj we define K = ^t] then in Bardeen coordinates K = p'^Cj, and 
{K, n) can be used in place of (w, p) as structure functions of /C flow. We have thus 
reduced the line element 12.11 of K, flow to 

(2.4) As^ =(l- ^1^^\ dt^ + 2K{p, z)dM</> - p'd<f' - e^^^"'^) {dp" + dz^) . 

We have seen, that in space-times given by the line element, world- lines of locally 
non-rotating observers are geodesies. 

At a given fixed point one can diagonalize the line element by a local linear 
transformation (local boost) of the base forms, say, di^ — > d0 -I- wdt. At the fixed 
point the line element reduces simply to dt^ — p^d(j)^ — e^^ {dp^ + dz^). The new 
base overlaps with the inertial frame of an observer momentarily co-moving with 
the locally non- rotating observer at that point. 

2.2. Equations for structure functions of K. flow. We have shown that for 
K, flow general metric 12.11 can be reduced to 12.41 In what follows we shall derive 
equations for K and p. Let E = G kT with T = VU ®U , U = n. Einstein's 
equation E = Q imply from 

E p-^\p —^ p. A, E^^-^yp ^ + M„ 



(EPp = -E\, EP, = E^'p) that 

(2.5) ^i,p=^p''^— — ^, ^1,^^-p-' 
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For C^ solutions the Schwarz identity p^pz = p.zp imposes on w the hnear elhptic 
constraint 

(2.6) Cj,pp + 3p-^LO,p + co,zz^0 ^ K^pp '^+K^zz = 0. 



As 



, e"^^ / K p 



equation E*' ^ = is satisfied identically. By calculating [i^p and /i^^ froni l2.^ and 
using W^ one may check the component E*■^ (then E'^^ — Kp^^E^^) reduces to 
E\ = -e"2^p2 ^(2,2^ j^ Qp._\^ - kP = 0, and finally, 

(2.7) KP--e-2V(<:^'p + ci'.). 

The other components of E vanish identically by symmetry. As energy density is 
negative definite /C flow is unphysical in the sense it cannot be made of ordinary 
matter. Once a solution of 12. 61 is found, which is simple task due to linearitv. 12.71 
gives the respective energy density and 12.51 can be easily integrated. Contribution 
to the total mass, from regions where ci) e C^, is negative and given for solutions 

by 

(2.8) Ms^ooth = -^ P^ {"^l + "^D dpdz 

anyway, we shall show that asymptotically flat space-times of IC flow must contain 
distributional sources of a nett positive mass —M smooth that cancels contributions 
(possibly infinite) from smooth components given by integrals as the above, as 
the asymptotic total mass, calculated as a surface integral at infinity, vanishes for 
asymptotically flat /C flows. 

Note, that elliptic equation 12.61 is the condition for critical points uj of energy 
functional 12.81 fwhen derived by integrating 12 . 71 with respect to the proper volume 
e^P^p of the surface of constant time, then Msmooth — J X'e^^'pdpd^dz) in any finite 
region of a hyper-surface of constant time bounded by a closed line in the plane 
[p, z) and such that uj is at least C^ (piecewise) in the region, namely 

SM[Oj]=0 ^ Cj,pp + 3-^ +Ci>^zz=0- 

For the solutions the maximum (minimum) principle applies. 

2.3. Some solutions. To obtain basic solutions we transform l2.6l to spherical coor- 
dinates p -^ rsin6', z — > rcosO. By substituting i^(r, 6*) = R{r)Y{9), the separation 
of variables gives R{r) — r~" or R{r) = r"+^, and the hypergeometric equation 
for Y{x), where x = cos 6. We assume here n G N and take only solutions that 
are analytic at a; = ±1, by which the multipole series is established. The general 



formula for z- antisymmetric and z-symmetric external multipoles that are solutions 
of 12. 61 and that give asymptotically flat space-times is 

(2.9) i.(p,z)H 0M--^^^^(™-^'-"^-'4^)' --1'2'3,.- 
(pi+J)^+U2 2F1 [m+^, -m; |, ^^J , m = 1,2,3,... 



For TO = one obtains K ~ ^ p^ + z^ and the monopole'* K = zj \J p^ + z^ which 
are non-flat asymptotically . The corresponding series of internal multipoles (sin- 
gular at infinity) is 



(2.10) K{p,z) 



(p2 + 22j™ ^p^ (^ _ 1 ^ _^. 1 ^ _^^ ^ m = 0, 1, 2, . . . 
z (p2 + 22^" ^p^ (^ _^ 1 ^ _^. 1^ -5^1^^ , m == 0, 1, 2, . . . 



2.3.1. An example. In the previous paper we found an example of asymptotically 
flat solution (note that equation for K function are identical but the models and 
line elements are essentially different) 

2p^z 



{z + a) J p^ + {z - a) + (z - a) J p"^ + {z + a) 



Ka{p,z) 



z — a+ \ p^ + iz — a) 
(2.11) . ^' - . V- V 




z + a+ J p^ + {z + a) 
which can be expanded in the series of external z-symmetric multipoles 
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Ka{p,z) ^ -^a ~- = a +— J- a 



where r = \/ p'^ + z^. Solution Ka has single extremum K{0, 0) — a and, apart 
from the singular points (0, a) and (0, —a), is everywhere smooth. The two points 
are singular since, for example, dzK as a function of z, for p — has a jump at 
z = ±a, a > 0. The corresponding angular velocity cja ~ P^^Ka has additional 
singularities, namely it is analytical everywhere with the exception of the closed 
segment p = 0, z G [—a, a]. The resulting space-time is asymptotically flat 

T^ f n\ 2 a^ . 2 ^ ~ t n\ 2 a^ a" sin"* 9 

Ka{r,e) ^ -—sm 9, uja{r,9) r^ -— Pair, 9) -— — , r ^ 00, 

6 r 6 r'^ 4r"^ 

and 

gtt^l- tfi sin^ 9, gt^ - ii^ sin^ (l + |^2^ (3 -H 5 cos 20) 

5^^ = -r2 sin2 9, g„^a^^-l-^ sin^ 0. 

Comparison with asymptotical expansion of the Kerr metric allows to determine 
total mass of the space-time, which is M = 0, and total angular momentum, which 
is J — a^/3. The result proves that the singular segment of solution uJa, located 
on the z-axis at —a < z < a, are sources of positive mass that cancel negative 
mass of the remaining regions {V < 0). In general, total mass of asymptotically 



We assume the same nomenclature for the ii" as in 2j where Euclidean mass density 
p~^ {K'^ + K\^ for the monopole of the van Stockum-Bonnor flow was spherically symmetric, 
even though, as one can check, the Euclidean mass density — p"*"^ ifi^o +'^i) f°'' ^'^'^ monopole 
u) = K/ p^ of the K, flow is not spherically symmetric. 



flat space-time can be calculated as surface integral over sphere S^ at infinity, and 
for the line element 12.41 

(2.12) M = lim — / •d^ = - lim — / / (r^ sin^ dCjdrUi) r^ sinMe* A d(p, 

r^oo 87r J r^oo StT J J 

where ^ = £,fj,dx^ and • is the Hodge operator. For r finite and sufficiently large 
the mass function is positive and behaves as M{r) ~ 4a'*/ (9r^) > 0. The function 
tends to as r — > 00, which confirms our statement. The singular segment is 
the space-time singularity as curvature scalar is proportional to T) which does not 
exist as a function on the segment, or more precisely, P is a distribution. Total 
angular momentum of an asymptotically flat space-time can be calculated as a 
surface integral over the sphere S^ at infinity, and for the line element 12.41 

(2.13) J ^ - lim -— -kdr] ^ - \im -— [r^ sin^ OdrCo) r^ smOde A d0, 

r— >oo IGtT J r^oc IGn J J 

and we again obtain J = a^ /3 for solution 12.111 As we have mentioned earlier, the 
specific angular momentum for any regular element of the fiow vanishes identically 
as trajectories of the fiow are world-lines of locally non-rotating observers. Another 
way of seeing this is to calculate the form d^drj, which measures angular momentum 
density when one take di = 0, and which as a general form vanishes identically for 
Cj satisfying eauation l2.6l In the singular regions where w is a distribution and does 
not satisfy 12.61 the contribution from d * drj to total angular momentum can be 
calculated from the above surface integral. 

2.4. Asymptotically flat space-times of /C flow are massless, some of them 
may have non-vanishing total angular momentum. The effect of vanishing 
of total mass should hold for all asymptotically fiat space-times of K. flow. Asymp- 
totical flatness for line element 12.41 means that w falls off like r~^ or quicker in 
spherical coordinates. If w ~ r~^ then such space-time may have nonzero total 
angular momentum. Vanishing of total mass is due to the presence of singularities 
that contribute positive mass that cancel contribution from regular regions of the 
space-times where energy density is negative definite. Similarly, the existence of 
asymptotically flat K, flows implies that the singularities they contain must have 
finite angular momentum equal to the total angular momentum of the space-times 
as the flow contributes no angular momentum in the regular regions. 

3. Conclusions 

We examined cylindrically symmetric and stationary dust flow along world-lines 
of locally non-rotating observers of which space trajectories are rings about common 
axis of rotation. Such flow can be constructed globally with the exception of regions 
where curvature singularities reside and that are additional sources of gravitational 
field. The fiow is differential (the shear tensor is non-zero), non-expanding (the 
dilation scalar vanishes identically) and locally does not rotate (the vorticity tensor 
vanishes identically). Local mass density of the flow is necessarily negative-deflnite 
which makes the flow nonphysical as far as usual forms of matter are concerned. 

There exist asymptotically flat space-times of the flow. Components of the flow 
move differentially on circular orbits with respect to an asymptotic stationary ob- 
server. The space-times have vanishing total mass and contain internal singularities 



where distributional sources of positive mass and of angular momentum are located, 
and this should be understood in the following sense. It turns out that a surface 
integral over the sphere at infinity and which reproduces total mass (of the Kerr 
black hole, for example), vanishes for the solutions. For radii sufficiently large the 
resulting mass function is positive and attains zero monotonically in the limit of 
infinite radius. As local energy density is negative definite in the regions where 
space-time is smooth, the internal singularities of curvature the solutions contain, 
must be distributional sources of positive mass (maybe infinite) of which contribu- 
tion to the total mass is screened by the regular regions, such the mass function 
is zero at infinity. This phenomenon is quite analogous with the screening of sin- 
gularities of negative mass by regions of positive energy density of asymptotically 
flat van Stockum-Bonnor flow. By construction of /C flow the specific angular mo- 
mentum per particle is zero, nevertheless, total angular momentum of a class of the 
asymptotically flat space-times may be non-zero and we gave an example in 12.111 
As the /C flow has no specific angular momentum, the total angular momentum, 
which is given by a surface integral over the sphere at infinity and which reproduces 
total angular momentum (of the Kerr black hole, for example), must be located in 
singularities of the space-times. 

Among other solutions, the model contains an infinite sequence of smooth asymp- 
totically flat multipolar solutions to which a class of other asymptotically flat solu- 
tions can be decomposed. There also exist an infinite sequence of the corresponding 
internal multipolar solutions that are not asymptotically fiat. The solutions are 
given in general by formulas 12.91 and 12.101 
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